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PREFACE 
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Shiffman. 
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R.  Covirant 


THEORY  0?  UiiJiihwATEH  KXPLOSIOH  BUBBLES 


Introduction 

Tne  theory  of  the  behayior  of  a  gas  bubble  produced 
b7  an  'jjidenTater  explosion  has  reached  such  a  3"a~e  of 
cornileteness  that  it  seens  desirable  to  give  a  unified 
acco'jnt  of  some  of  the  nain  results,   'rfhen  an  explosive 
is  detonated  ^jnder  water,  a  shock  Tave  is  first  enitted 
and  then  the  gaseous  products  of  za.e   explosion  expand 
Tinder  the  influence  of  their  high  internal  presa'Jire.   Be- 
cause of  its  inertia,  the  gas  bubble  overexpands  to  a 
verj  low  press'ore  and  then  zae   hydrostatic  press-J-re  of 
the  water  recoaprssses  the  bubble  to  nearly  its  original 
size.  At  this  stage  the  bubble  starts  expanding  again 
and  a  press"CLre  pulse  is  eni-ted.   This  process  of  expan- 
sion and  contraction  may  occur  five  or  six  tines  before 
the  bubble  breaks  up  an^^  dissolves  or  escapes  fron  the 
water. 

I^ne   behavior  of  ~he  bubble  is  affected  by  the  pres- 
ence of  3'Jurfaces,  such  as  the  cotton  of  the  water,  the 
walls  of  a  tan>,  the  target  or  the  air-water  s'urface. 
!;cm=ll7,  the  bubble  would  tend  to  nove  upward  since  it 
is  lighter  than  the  surrounding  water.   However,  the  free 
surface  of  the  water  repels  the  bubble  while  a  rigid  s^-lt- 
face  of  any  kind  such  as  a  wall  or  a  'zotzaa.   attracts  it. 
In  sone  cases  these  effects  nay  be  so  strong  that  the 
bubble  actually  noves  downward. 

Besides  influencing  the  notion  of  the  b-ubble,  the 
surfaces  also  change  the  period  of  oscillation  snl  nodify 
'*  the  peak  press'ores  produced.   The  theory 


presented  here  explains  these  effects  and  gives  a  quanti- 
tative estimate  of  their  size.  The  vsQ.ues  predicted  for 
the  period  show  excellent  agreement  with  those  obtained 
experimentally,  while  the  values  predicted  for  the  peak 
pressure  and  the  distance  moved  by  the  bubble  show  only 
fair  agreement. 

An  important  application  of  the  theory  is  the  deter- 
mination of  the  amount  of  energy  left  in  the  bubble  after 
the  shock  wave  has  passed.  Also  the  energy  in  successive 
bubble  oscillations  can  be  foxrnd.   There  is  a  large  amo\int 
of  energy  dissipated  in  the  transition  from  shock  wave 
stage  to  bubble  motion  and  also  between  successive  bubble 
oscillations,  which  cannot  be  explained  on  the  basis  of 
the  energy  radiated  by  the  pressure  pulse.   The  explana- 
tion of  this  large  energy  dissipation  is  still  unknown. 

Outline 

As  developed  in  this  paper  the  theory  is  an  exten- 
sion of  that  given  in  AMP  Report  37. IR,  studies  on  the 
Gas  Bubble  Resvilting  from  Underwater  Explosions;   On  the 
Best  Location  of  a  Mine  Near  the  Sea  Bed.   There,  the 
motion  of  the  bubble  in  the  presence  of  a  rigid  bottom 
was  Investigated  and  it  was  shown  that  the  exact  theory 
could  be  successfully  approximated  by  the  addition  of  a 
term  to  the  kinetic  energy.   In  this  paper  we  show  that 
the  effects  of  surface,  bottom,  walls,  targets,  etc.,  can 
all  be  approximated  in  the  same  way  by  the  addition  of  a 
suitable  term  to  the  kinetic  energy.   The  evaluation  of 
this  term  depends  upon  the  solution  of  an  "electrostatic 
problem."  In  Section  IV  we  work  out  in  detail  the  case 
of  a  bubble  between  a  free  siorface  and  a  bottom.   Other 
cases  can  be  treated  in  the  same  way. 

Section  I  presents  a  collection  of  formulas  and  a 
stunmary  of  methods  which  can  be  used  to  predict  the  pe- 
riod of  oscillation  of  the  bubble,  the  distance  its 


center  moves  during  the  first  oscillation,  the  maximum 
and  minitnvun  radius  of  the  bubble  and,  finally,  the  peak 
pressure  emitted  by  the  bubble.   The  formulas  are  given 
in  terms  of  certain  integrals  which  can  be  evaluated  by 
the  method  given  in  Report  37. IR.  Nvunerical  integration 
of  the  differential  equations  of  motion  of  the  bubble  is 
conipletely  lonnecessary.  For  convenience  Figures  1-6  con- 
tain graphs  of  these  integrals  for  the  most  frequently 
occurring  values  of  the  parameters. 

Section  III  contains  a  discussion  of  these  formulas 
and  a  short  indication  of  their  proof.  A  careful  analy- 
sis is  made  of  the  dependence  of  the  parameters  in  the 
bubble  motion  upon  the  properties  of  the  explosive.   It 
is  shown  that,  by  a  study  of  the  periods  of  bubbles 
placed  at  different  depths,  it  is  possible  to  determine 
the  amount  of  energy  left  in  the  bubble  after  the  shock 
wave  has  passed  and  also  to  determine  the  exponent  in  the 
equation  for  the  adiabatic  expansion  of  the  explosive. 
This  seems  to  be  one  of  the  very  few  methods  by  which 
this  exponent  can  be  found. 

A  similar  procedure  can  be  used  to  determine  the 
amount  of  energy  left  in  the  second  bubble  oscillation. 
The  experiments  indicate  that  only  about  16  percent  of 
the  original  energy  of  the  explosive  remains.   CaJ-cula- 
tions  show  that  the  energy  radiated  by  the  pressure  pulse 
emitted  by  the  bubble  at  minimum  size  is  not  large  enough 
to  explain  the  energy  loss.  The  explanation  of  this  high 
dissipation  of  energy  is  one  of  the  major  imsolved  prob- 
lems of  the  theory. 

As  was  mentioned  before.  Section  IV  contains  a  solu- 
tion of  the  "electrostatic  problem"  equivalent  to  the 
problem  of  a  bubble  placed  between  a  rigid  bottom  and  a 
free  surface. 

In  Section  II  the  theory  is  applied  to  the  analysis 
of  some  experimental  data  obtained  at  Woods  Hole  by  Arons 


and  his  co-workers.   Three  hxindred  grams  of  Tetryl  were 
detonated  at  varying  depths  below  the  surface  in  water 
23.5  feet  deep.   The  period  of  oscillation,  the  peak 
pressure  and  the  distance  the  bubble  moved  were  measured. 
The  agreement  between  theory  and  experiment  as  regards 
periods  is  excellent;  as  regards  pressure  and  distance 
the  bubble  moves,  the  agreement  is  only  fair. 

Since  the  agreement  as  regards  periods  is  one  of  the 
outstanding  successes  of  the  theory,  it  seems  worthwhile 
elaborating  on  it.  Let  T  be  the  period  of  oscillation  of 
the  bubble  at  a  depth  H  feet  beneath  the  surface.   In 
Section  II  it  is  shown  that  T(H  +  33)  '      varies  linearly 
with  ^J-.  (H  +  33)"  '      where  9?-,  is  a  complicated  fxonction  of 
H,  so  that  we  have 

T(H  +  33)^/^  =  a  +  (3(?)^(H  +  SS)"-"-/^ 

where  a  and  p  are  constants  independent  of  H.   It  is  also 
shown  that  a  =  C-j^E  ^"^  and  p  =  ^^       where  E  is  the 
amount  of  energy  left  in  the  bubble  after  the  shock  wave 
has  passed  and  C-,  and  Cg  ar©  two  constants  depending  upon 
the  exponent  in  the  equation  for  the  adiabatlc  expansion 
of  the  explosive. 

Prom  the  experimental  data,  T(H  +  33)^   and 

-1/3 
^•j^(H  +  33)   '   are  calculated  and  then  the  constants  oc   and 

p  are  determined  by  the  method  of  least  squares.   Prom  the 
values  of  a  and  p  the  amount  of  energy  left  in  the  bubble 
and  the  adiabatic  constant  can  be  determined.   In  the  par- 
ticular case  considered  in  Section  II  it  is  found  that 
about  48  percent  of  the  original  energy  of  the  explosive 
is  left  in  the  bubble. 

Assumptions 

As  in  Report  37. IR  we  shall  idealize  the  problem  by 
making  the  following  assumptions:  . 


1.  The  water  is  an  ideal  incompressible  fluid. 

2.  The  bubble  remains  spherical  in  shape. 

3.  The  gas  inside  the  bubble  expands  adiabatically. 
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Section  I 
Summary  of  Formulas 


Units  of  Length  and  Time 

Let  A  be  the  radius  of  the  bubble  at  time  T,  B  the 
distance  of  Its  center  from  the  bottom,  and  H  the  dis- 
tance from  the  surface.  We  put  Z  =  H  +  Z",  where  Z"  Is 
the  height  of  the  water  colxmn  whose  pressure  equals  at- 
mospheric pressure,  so  that  the  hydrostatic  pressure  at 


h 


B 
1 


UatzTSur/cii 


ce 


bofiom 


the  center  of  the  bubble  is  pgZ. 

Let  E  be  the  amount  of  energy  left  in  the  bubble 
after  the  shock  wave  has  passed  and  let  M  be  the  mass  of 
the  explosive.  We  write  E  =  rQM  so  that  rQ  Is  the  amoxint 
of  energy  per  xmit  mass  of  explosive  left  in  the  bubble 
after  the  shock  wave  has  passed. 

All  the  formulas  will  be  given  in  dimensionless 
terms.   The  unit  of  length  is 


(1.1) 


L  =  (5E/4TrpgZ^)^/^ 


where  Z^  is  the  value  of  Z  at  the  time  of  detonation. 
The  vmit  of  time  is 


(1.2)  C  =  L(2gZQ/3)"^/^. 

Note  that  the  xmlt   of  velocity  is  proportional  to  the 
square  root  of  Z   »   for 

(1.3)  L/C  =  (2gZy3)^/2. 

Parameters  in  the  Bubble  Motion 

If  the  pressure  and  volume  of  a  unit  mass  of  the  gas 
formed  by  the  explosive  are  connected  by  the  adiabatic 
relation: 

(1.4)  PV"''  =  K,  ' 

then  the  internal  energy  of  the  bubble  when  expressed  in 
non-dimensional  terms  will  depend  upon  a  parameter  k 
which  is  defined  as  follows: 

(1.5)  k  =  K(E/M)"'^(pgZQ)^"V(Y  -  1). 

Since  y  -  1  Is  very  close  to  zero,  the  value  of  k  will 
hardly  change  when  Z  varies  slightly.   In  most  cases, 
therefore,  we  shall  be  able  to  treat  k  as  constant  inde- 
pendent of  Z  . 

The  effects  of  free  surface,  gravity  and  bottom  are 
expressed  with  the  help  of  a  function  §0-,  .  We  put 

(1.6)  p^=   -[f (x)  +  log  2]/(H  +  B) 

=  -^(x)/H 

where 

(1.7)  ^=1-^ 

and 

00      ,  vn 

(1.8)  f  (x)  =  2x  z: ^-^ K 

o  (2n  +  1)'^  -  -xr 
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00      _/  ^n 


fr.„.  -    X     .  ,,     ^^r^ ^i-r .logs 


and 


(1.9)  «:/{x)  =  ,p^  +  (1  -  X)[Z:  ^^ 5  +  ^^^1 

Tables  and  graphs  of  the  functions  f{x)  and  t^ (x)  are 
given  in  the  appendix  (see  Figiires  13,  14  and  15  and 
Tables  3  and  4) .  Note  that  if  the  water  is  infinitely 
deep,  B  =  00  so  that  from  (1.7),  (1.9)  and  (1.5)  we  have 

(1.10)  X  =  1,   <^(x)  =  1/2,  <p^   =  -1/2H. 

Bubble  Quantities  at  Time  of  Maximum  Size 

Small  letters,  a,  b,  h,  etc.,  will  be  used  to  repre- 
sent the  non-dimensional  values  of  the  quantities  repre- 
sented by  the  corresponding  capital  letters.   The  sub- 
script zero  attached  to  these  quantities  will  indicate 
their  value  at  the  beginning  of  the  bubble  motion,  the 
subscript  one  their  value  at  the  time  of  maximum  size, 
and  a  bar  attached  to  them  will  Indicate  their  value  at 
the  time  of  recompression  when  the  bubble  size  is  a  mini- 
mum. 

The  non-dimensional  time  from  the  beginning  of  the 
bubble  motion  to  the  stage  of  maximum  size  is 

(1.11)  *1  =  ^1  +  ^  l2^^1' 
where 

(1.12)       i^  =  r^      ^  ^^ 


a. 


^1  -  a»  -  ta-'(v-lV 


and 


(1.13)        ig  =  r^        a^/^  ^^r— . 

'  i         ^1  -  a^  -  lca-^^^-1)' 
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Here,  a  and  a,  denote  the  smallest  and  largest  roots, 
respectively,  of  the  denominator.   Graphs  of  I,  and  Ig 
for  the  most  frequently  occurring  values  of  k  and  y  are 
given  in  Figures  1  and  2. 

The  vertical  momentum  at  the  time  of  maximum  size  is 

(1.14)  s^  =  L  -^(13  -  I  Lp^l^)    +  2^^(15  +  ^  L^^^Ig) 
where 

(1.15)  I3  =  r^   a^/^^l  -  a^  -  ka-5(v-i)'  da, 

(1.16)  I,  =  /"I  B^/^il  -   a=  -  ^^-■■"■■^-'■y   da. 


■4 

^o 


^^1       a^/2  ^^ 


^""^       '^-(  irrrr^m^' 


(1.18)       I 


^1      a"/2  da 


'       i       \/l  -  a'  -  to-'(v-^)'  ■ 


The  maximum  radius  of  the  bubble,  a-,,  is  the  root 
near  1  of  the  equation: 

(1.19)  *    a^  +  ka'^^*^-^^  +  3s^/2a^  =  1. 

A  good  approximation  to  a-,  is  given  by  the  formula: 

(1.20)  a^  =  1  -  I  -  (Y-jJ:  -.  I)k2  -  I  s2. 

Bubble  Qtiantities  at  Time  of  Minimum  Size 

The  time,  t,  from  the  beginning  of  the  bubble  motion 
to  the  stage  of  minimiom  size  is  twice  the  time  to  the 
stage  of  maximum  size,  that  is, 

(1.21)  t  =  2t^. 
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The  vertical  momentum,  i,  at  the  time  of  minimum 
size  is  twice  the  vertical  raoraentum  at  the  time  of  maxi- 
miim  size,  that  is, 

(1.22)  s  =  2s^. 

The  minimtim  radius  of  the  bubble,  a,  is  the  root 
near  zero  of  the  equation: 

(1.23)  a^  +  ka"^^"^'^^  +  3i^/2a^  =  1. 


If  we  put 

(1.24) 

5  =  j,l/(Y-l),        , 

=  -bl^/2% 

and 

/ 

(1.25) 

i^  =  ^u. 

then  u  is  the  root  of  the  equation: 

(1.26)  u  =  u^"*''  +  £. 

A  graph  of  u  for  various  values  of  e  and  y  is  given  in 
Figure  7. 

The  non-dimensional  distance  the  bubble  moves  down- 
ward is  given  approximately  by  the  formula: 

(1.27)  Ab  =  3i(I^   +  \  LjPj^Ig) 
where 

a-^/2  da 


(1.28)         I„   =     /     ■■• 

^        i         ^1   -  a^  -  ica-^^V-l)    -   3s-2/2a^ 

and 

(1.29)   iq  =  r^ .  ^"^^^  ^^  =,  . 

^       i         \/l  -  a^  -  ka-^^-^-l)    -  3i2/2a=^' 
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Tlie  non-dimensional  peak  pressiore  at  the  time  of  mln- 
imvm  size  is 

(1.30)  (a^a)-   =  r^/^  u-5/3[£2  ^  Zj^  ^a-Y^^i  ^  ^^yl 
where 

(1.31)  u  =  a^  5"-'-. 

Actual  Values  of  Bubble  Quantities 

If  f  is  the  actual  period  from  the  beginning  of  the 
bubble  motion  to  the  time  of  wiiniTtrnm  size>  then 

T  =  CB 

or  using  (1.21),  (1.11),  (1.1)  and  (1.2)  we  find  that 

(1.32)  T  =ctz:^/^  +  ^f.K^^^ 
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where 


(1.33)  a  =  2l3^(2g/3)"^/^  (2Trp)"^/^  E^/^, 

(1.34)  p  =  l2(2g/3)-'^/^  (2ttp)-2/^  E^/^. 
Note  that 

(1.35)  TZ^/^  =  0.  +  P^iZ;^/^, 

so  that  fz^''   is  a  linear  fvuiction  of  <P-,ZZ  ^    >      In  a  later 
O  X  o 

section  we  shall  use  this  formula  to  find  a  and  p  and  then 
the  value  of  E. 

The  maximxmi  and  minimum  radii  are  fovmd  from  formulas 
(1.20)  and  (1.23)  after  multiplication  by  the  tinit  of 
length,  L,  given  in  (1.1). 

The  peaJc  presstire  of  the  bubble  at  a  distance  R  from 
the  center  of  the  bubble  is 


12 


(1.36)  P  =§  psZ^.^(a^a)- 

where    the  value   of    (a^a) *    is   given  by  formula    (1.30). 

At   the   time   of  majcimim  size   the  pressure   in  the  bub- 
ble  is 

P  3(1   -   5aJ) 

(1.37)  Pi=|pSZo  —2 
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Section  II 
Analysis  of  Experimental  Data 


In  this  section  some  experimental  data  obtained  by 
A.  Arons  and  his  co-workers  at  Vi/oods  Hole  Oceanographic 
Institution  will  be  analyzed.   It  will  be  shown  that  the 
theoretical  values  of  the  period  agree  very  well  with  the 
experimental  data  while  the  values  for  pressure  and  mi- 
gration do  not  agree  as  well. 

The  experiments  consisted  of  exploding  three  hundreo 
grams  Tetryl  at  various  depths  beneath  the  surface  in 
water  23,5  feet  deep.   Table  1  shows  the  values  observed 
for  the  period  of  bubble  oscillation. 

Before  these  data  can  be  compared  with  the  theoretical 
results,  the  value  of  E  must  be  known.  We  shall  find  E 
by  fitting  a  straight  line  to  the  points  obtained  when 

fz  '     is  plotted  against  (p-,Z'  '     .   Table  2  gives  the  val- 

-  5/6  -1/3 

ues  of  TZ    and  the  corresponding  values  of  9-{L~        . 

A  straight  line  is  fitted  to  these  values  by  the 

method  of  least  squares.   The  equation  obtained  is 

(2.1)  TZq/^  =  3.876  -  16.96  <P-^~^^^ ' 

pigvire  8  shows  how  closely  the  straight  line  fits 
the  data.   This  closeness  of  fit  is  a  confirmation  of  the 
theory. 

When  equation  (2.1)  is  compared  with  (1.35),  we  find 
thet 

(2.2)  a  =  3.876,    p  =  16.96. 
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Table  2 


H 



^.<'' 

H 

TZ^/« 

H-l''' 

H 

r^r 

f.^-o''' 

3.5 

3.158 
3.160 

.0428 

7.0 

5.524 
3.554 

.0198 

15.4 

3.830 

.000620 

3.170 

3.560 

16.4 

3.893 

-.00210 

4.0 

3.223 

.0372 

3.571 

17.3 

3.969 

-.00525 

3.239 

8.0 

3.575 

.0166 

18.5 

3.938 

-.0108 

3.247 

3.617 

3.263 

3.621 
3.634 

18.8 

3.971 

-.0122 

4.5 

3.290 

.0332 

19.0 

3.915 

-.0140 

3.296 

10.0 

3.680 

.0121 

3 .  933 

3.307 

3.634 
3.696 

19.1 

4.032 

-.0148 

5.0 

3.353 

.0292 

3.707 
3.710 

3.357 

19.2 

3.930 

-.0157 

3.363 

12.0 

3.722 

.00710 

19.6 

3 .  925 

-.0138 

5.5 

3.419 
3.430 
3.438 

.0262 

3.732 
3.735 
3.746 

19.7 

3.963 
3.996 
4.167 

-.0197 

6.0 

3.442 

.0238 

3.796 

19.8 

3.917 

-.0204 

3.448 

14.0 

3.788 

.00385 

3.953 

3.452 

3.795 

3.969 

6.5 

3.516 
3.527 
3.531 

.0216 

3.817 

20.2 
20.4 

4.091 
4.174 

-.0264 
-.0288 
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Formulas  (1.35)  and  (1.34)  shov/  that  a  and  6  depend 
on  the  value  ol'  E  and  of  I^  and  Ig.   If  we  consider  the 
ratio  ct  /p  we  have 

2  I  2 

(2.3)  |-  =  4l2i-l(2g/3)-^/2   or   ^  =  1.155  ^ 

p 
which  is  independent  of  E.  A  graph  of  Ij^/lg  jfor  different 

values  of  k  and  y  is  given  in  Figure  9. 

Three  physical  quantities  are  needed  to  fully  describe 
the  bubble  motion — the  values  of  K,  Y  and  E.   Since  we  have 
only  two  constants  a.  and  p  determined,  we  must  assume  the 
value  of  one  of  the  three  quantities.   We  shall  assume  that 
Y  =  1.25,  the  value  proposed  by  Jones  for  T.N.T.,  and  then 
by  the  use  of  Pigiire  9  we  find  that  k  =  .23.   Prom  this 
valxie  of  k,  using  Figiires  1  and  2  and  formulas  (1.35)  and 
(1.34),  we  find  that  rQ  =  500  calories/gram.   Since  the 
detonation  energy  of  Tetryl  is  1060  calories/gram,  this    ' 
value  for  rQ  is  another  confirmation  of  the  theory. 

Since  we  know  the  value  of  E  =  rQI/I,  we  can  use  (1.1) 
and  get 

(2.4)  L  =  15.88  (w/Zq)"""/^ 

where  L  is   in  feet,   W  in  poimds   and  Z      in  feet,    and 

C   =  2.997  \?-^\''^^^ , 
o 

v-1 
Since  k  varies  with  depth  as  Z    ,  the  value,  k  =  .23, 

is  really  an  average  value  over  the  range  Z  =  33  to 
Z  =  56.5.   If  we  assume  that  this  value  of  k  corresponds 
to  a  depth  halfv/ay  between  the  bottom  and  the  surface  of 
the  water,  we  find  that  k  varies  from  .215  to  .244.   This 
change  in  k  can  be  neglected  in  the  calculation  of  the 
period,  momentum  and  migration.  However,  in  the  calcula- 
tion of  the  peak  pressiore  it  must  be  taken  into  account 
since  the  pressvire  varies  as  k~  /'*''"' . 
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Pigiire  iu  shov;s  how  the  periods  calculated  rrom  the 
equation  (2.1)  compare  with  the  observed  results.   Tne 
agreement  is  excellent. 

Using  k  =  .23  and  formulas  (1.27)  and  (1.28)  the  mi- 
gration v/as  found  by  numerical  integration  of  I_  for  each 
value  Z  .   The  calculatea  results  are  compared  v/ith  the 
observed  values  in  Figure  11.   The  agreement  seems  to  be 
fair. 

The  peak  press\ire  is  calculated  froi;i  formultts  (1.56) 
and  (1.50).   If  vie   use  the  values  of  k  and  y  aetermined 
by  the  period  measurements,  the  calculated  pressures  are 
about  one-half  the  observea  pressure.   Theoretically,  it 
would  be  possible  to  use  one  experimentally  determined 
value  of  pressure  with  the  tv/o  constants  a  and  p  to  solve 
for  the  three  quantities  rQ,  k  and  y«   This  was  done  but 
it  did  not  lead  to  consistent  results.   It  seems  likely 
that  the  value  of  y  at  the  time  of  minimum  size  may  be 
different  from  the  value  of  y  during  the  time  of  expan- 
sion and  contraction. 

Despite  this  difficulty  with  the  magnitude  of  the 
pressvire,  the  variation  of  pressure  with  depth  is  given 
approximately  by  formula  (1.50).   If  we  multiply  the  cal- 
culated pressxxres  by  a  constant  chosen  so  that  the  calcu- 
lated pressure  agrees  with  that  experimentally  observed 
at  the  depth  of  18.5  feet,  vie   get  Figure  12. 
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Section  III 


Derivation  or  formulas 


Consiaer  the  motion  of  a  gas  bubble  in  water  of  fi- 
nite depth,  talking  into  accovmt  the  effects  of  the  bottom, 
the  free  sxorface  and  gravity.  As  before,  we  let  A  be  the 
radius  of  the  bubble  at  the  time  T,  B  the  distance  of  its 
center  from  the  bottom,  and  H  the  distance  from  the  sur- 
face,  we  put  Z  =  H  +  z"  where  z""  is  the  height  of  the 


H 


B 

1 


water  column  whose  pressiire  equals  atmospheric  pressure, 
so  that  the  pressure  at  the  center  of  the  bubble  is  pgZ . 
Just  as  in  Report  57. IR  we  find  that  if  we  use  coor- 
dinates R  and  0  to  describe  the  motion  of  the  water,  the 
velocity  potential  describing  the  flow  is 


{-6.1) 


.^aV.aVoo^,,^ 


v/here  $-,  is  the  "image"  potential  necessary  to  satisfy 
the  bounaary  conditions  on  the  bottom  and  on  the  surface, 
The  primes  indicate  time  derivatives. 

By  classical  hydrodynamics  we  can  show  (see  Report 
37. IR)  that  the  kinetic  energy  of  the  water  is  given  by 
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(3.2)  2^pa^^La'''^(1  +  P.p^)    -   2a^^2^'b'  +  ^  b'''^(1  +  k^<p^)] 

where  cp   ,   f^   ana  cp.     are  functions  of  A»  B  and  H  which 
represent  the  effect  of  the  surface  and  the  bottom.   Note 
that  9P,  ,  <p^,   9?„  are  of  the  dli-nonsions  (length)"", 

(length)"  ,  (length)"'  respectively. 

Ihe  potential  energy  of  the  displaced  water  is  equal 
to  the  volume  of  the  bubble  multiplied  by  the  hydrostatic 
pressure  at  the  center  of  the  bubble,  that  is,  to 
4TrA^pgZ/3 . 

Assume  the  pressxire  and  volxime  of  one  unit  mass  of 
the  gas  in  the  bubble  under  an  adiabatic  change  are  con- 
nected by  the  relation 

(3.3)  PV*^  =  K. 

If  M  is  the  mass  of  the  explosive,  then  for  the  actual 
bubble 

PV*^  =  KI^*^ 

and  the  internal  energy  of  the  gas  when  the  bubble  has 
radius  A  is 

G(A)  =  J  Pdv  =  iai'^v"'-"V(r  -  1) 


(3.4) 


YJ^^d  iTA^)^"V(r  -  1) 


By  our  assumption  the  total  energy,  E,  left  in  the 
bubble  after  the  shock  wave  has  passed  is  equal  to  the 
sum  of  the  kinetic  energy  of  the  water  and  the  potential 
energy  of  the  displaced  water  and  the  internal  energy  of 
the  gas .   V/e  have ,  then 

2ttpA^[A'^(1  +  Afi)  -  2a^^2^'b'  +  \   b'^(1  +  A^^^)] 

(5.5)  .         , 

+  ^  TTA  pgZ  +  G(A)  =  E. 
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If  we  use  L  and  C  defined  in  (1.1)  and  (1.2)  as 
vmits,  equation  (3.5)  can  be  written  in  terms  of  non- 
dimensional  quantities  as  follows: 

a^[a^(l  +  aLjP3_)  -  2.Q^^Z?<f>^   +  ^  b^(l  +  a^^so^)] 

+  a'^zz^-^  +  ka  ^^^  ^^    =  1, 

v^here  ka"^^*^"'''^  is  the  internal  energy  of  the  gas  ex- 
pressed in  non-dimensional  quantities,   k  is  defined  in 
equation  (1.5) . 

It  was  shown  in  Report  37. IR  that  the  period  can  be 
found  by  assijraing  b  is  constant  so  that  b  =  0.  Equation 

(3.6)  reduces  to 

(3.7)  a^a^(l  +  aL^^)  +  a^  +  ka"^^*^""^^  =1. 

If  we  approximate  (1  +  aL^^)  '   by  1  +  aLf?-j^/2,  equation 

(3.7)  can  be  solved  for  a  and  we  get 

a^/^(l  +  aL9?T/2)da 

(3.8)  ,     ,      :,   ,,,  =  dt. 

V'l  -  a^  -  ka-^^-^-l) 

Let 


(3.9) 


(3.10) 


where  a  and  a,  are  respectively  the  smallest  and  largest 
zeros  of  the  denominator.  V/e  find  then  that  the  time 
from  beginning  to  maximum  is 

(3.11)  ^1  =  II  +  ^  l2^^1' 


I,  =  r^ 

a'/2  ^^ 

lo  =  T' 

a5/2  da 

2   J    X- 

.  a'  -  ka-'C^-l)' 
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which  Is  formula  (1.11).   The  total  period  Trom  beginning 
to  minimum  will  be  twice  t-,  . 

Momentum  of  the  Bubble 

The  vertical  momentum,  s,  of  the  bubble  is  equal  to 
a  "h/Z.      Using  the  Lagrangian  equation  associated  with  the 
b  coordinate,  we  find  that 

(3.12)       ^(aV3)  =  a^^^L  ^  +  a^z^^ 
Integration  gives  the  value  of  s  at  the  maximumi  as 

s^t  =  2[f  ^   a^a^L  -^  dt  +  z^'^  J    ^   a^dt]  . 

o  o 

Using  equations  (3.7)  and  (3.8)  and  the  fact  that 

adt  =  da 
we  express  s,  as  integrals  over  a.   We  have 


s,  =  L 


1   "  db 

a 


f^   a6/2(l+aL^^)-l/2l/l.a5-ta-2(Y-l)-  ^^ 


,   na,   a^'^d  +  aLa>-,/2)da 
(3.13)      +  z„   /  ■  i 

t'l  -  a'  -  ka-SC<-l)' 


o 
a 


o 


where 


22 


I^  =  jT'l  a^/^il   .   a^  -  ka-^(^-^)'  da, 


^o 


(3.14) 


I^  =  f^  s^/^h  -   a^  -  ka-^^^-1)"  da, 

ri       a^/^  da 
o 

r^  a^^/g  da 

'   i    /l  -  a^  -  ka-3C^-l)'  ■ 
o 

These  equations  are  the  same  as  (1.14)-(1.18) . 

The  momentxim  can  be  used  to  give  us  an  approximation 

to  the  energy  equation  (3.6)  which  takes  into  account  the 

fact  that  b  does  change,   we  shall  neglect  the  terms 

2  2         3  3 
a  L  jDp  and  a  L  ^^  which  are  of  second  and  third  order, 

respectively,  as  compared  to  the  first  order  term  aL^-,  . 

Now,  since 

b  =  3s/a^, 

equation  (3.6)  becomes 

(3.15)  a^a^(l  +  aL^^^)  +  a^  +  ka"^^"^""^^  +  3s^/2a^  =  1. 

If   s  were  knov/n  as   a  function  of  a,    this   could  be   inte- 
grated: 

p  a'^/^(l  +  aL5DT/2)da 

(3.16)  /     ■  ^        ^  .  =  t. 

^    h   -  a'>  -  ka-5(*^-l)    -   3s2/2a^ 

The  vertical  distance  the  bubble  moves  can  be  found 
from  this  formula  since 

S  =  3s/a^ 
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so  that 

p  -z  n        sa""^/^(l  +  aL^-,/2)da 

(3.17)   Ab  =  3  /  sdt/a"^  =  3  /  ,  ■ .   -.  x     ^   ^'  • 

^  ^  \/l-a^-lca-=^(*^-l)-5s2/2aS 

Since  s  is  given  as  the  indefinite  integral  of  the  right 
side  of  (3.12),  Ab  is  really  a  double  integral  but  we  can 
change  it  to  a  single  integral  by  the  follov/ing  argvunent: 
Dirring  the  first  part  of  the  bubble  motion,  that  is, 
until  the  bubble  reaches  its  maximum  size,  s  increases 
from  zero  to  s, .  During  the  second  part  of  the  motion, 
from  the  bubble  maximum  to  the  bubble  minimum,  s  increases 
from  S-,  to  s  =  2s^.   The  momentum  at  any  time  during  the 
second  half  of  the  motion  is  i  minus  the  momentum  devel- 
oped as  the  bubble  contracts  from  its  given  size  to  its 
minimiim  size.   Because  of  the  symmetry  of  the  bubble  mo- 
tion, this  minus  momentum  is  the  same  as  the  momentvim  de- 
veloped when  the  bubble  increases  from  its  initial  size 
to  the  given  size.  We  may  then  conclude  that  the  dis- 
tance moved  during  the  bubble  expansion  just  balances  the 
distance  which  would  be  produced  by  the  minus  momentum  so 
that  the  total  distance  moved  is  given  by  one  term 

,a      a"^/^(l  +  aL^  /2)da 
Ab  =  3s  ' 


(3.18)  ^1   ^1  -  a^  -  ..-=<-^)  -  3iW 


=  3s(I^  +  ^  L^i^e)' 

-3/2 


727Z^  ' 


-  3sV2a' 


1 

/^  a~~' ^  da 

,   Vl  -  a^  -  ka-^^^-^-l"-  3s-2/2a^ 


This  proves  (1.27) . 
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Radius  of  the  Bubble 

The  maximum  radius  of  the  bubble  is  found  by  putting 
a  =  0  in  equation  (3.15)  and  solving 

(3.21)  a^  +  ka'^^*^"^^  +  3s^/2a^  =  1 

for  the  root  near  one.   If  we  expand  a  in  powers  of  k  and 
solve  (3.21)  we  find  that 

3g2 

(3.22)  a^^  =  1  -  I  -  (^^-§-i  +  |)1^^  -  "2^  • 

This  is  the  same  as  formula  (1.20).  ' 

To  find  the  minimum  radius  of  the  bubble  we  put 
a  =  0  in  equation  (3.15)  and  consider  it  at  the  time  of 
minimum  size.   The  equation  reduces  to 

(3.23)  a^  +  ka"^^'^"-'-^  +  Ss^Sa^  =  1 

and  the  root  near  zero  is  the  minimvim  radius .  Note  that 

(3.23)  is  the  same  as  (3.21)  except  for  the  value  of  the 
momentum.  The  equation  simplifies  if  we  put 

(3.24)  ^  =  k^/^'^-^) ,    i.^   =tvL,         6  =  3iV23^ 

3 
and  neglect  that  a  term,  since  a  is  near  zero.  We  have 

U-('^-l)  +  £U-1  =  1 


or 

(3.25)  u  =  u^""^  +  6. 

Peak  Pressure 

The  pressure  in  the  water  is  given  by  Bernoulli's 
Law  as 


p   at  ^^  *^  • 
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The  second  term  can  "be  neglected  compared  to  the  first. 
Using  (3.1)  we  find  that 

P  _  (a^a')'  ,  (A^b')'  cos  e  ^  ... 

For  moderate  distances  from  the  bubble  only  the  first 
term  is  important,  so  that 

The  value  of  (a  a)  '  can  be  foiond  by  multiplying 
(3.15)  by  a  and  then  differentiating.  This  gives 

2a^a(a^a)*(l  +  aL^^)  +  a^a^aLf^  +  4a^a 

+  k(4  -  3Y)a^"^'^4  -  3s^a"^a  =  a. 

Since  the  peak  pressure  will  be  found  at  the  time  of  min- 
imum size 
result  is 


imum  size,  we  may  put  a  =  0  and  solve  for  (a  a)*.  The 


1  ^  3i^   k(4  -  5y)  -1-3y 

,  2-x-    2a'^ 2ar  '^ 

(a  a)   =  


I  + 


aL^-i 


(3.27)  o 

^^      +  5k(Y  -  1)  £l-3Y 


45^ 


1  + 


aLf-j_ 


if  we  make  use  of  equation  (3.23).   Using  the  notation  of 
(3.24)  equation  (3.27)  becomes 


^o2•^.  _  1  r3e  .  3(Y  -  1)  ,,1-Yt 
(a  a)   -  r^t^  +  5    ^   J 


(3.28) 


-  3(Y  -  1)  r  .    £  _  .   +  n^-Yn 

-  2^V3  ^2/3t(Y  -  l)u  ^  ^   J  • 
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Notice  that  since  u  Is  near  one  the  bracket  does  not 
change  much  as  k  changes . 

We  may  write  formula  (3.26)  for  the  pressure  in  a 
more  convenient  form  as  follows : 

c 


(3.29) 


^2/3  ^S/li  Wy  -  l)u  ""^        ^R 


by  formulas  (1.3)  and  (3.28). 

Replace  L  by  its  value  given  in  (1.1);  then 

(Y  -  DpgZ      ^       1-Yw   5E   ^1/5 

This  formula  can  be  combined  with  the  formulas  (1.33)  and 
(1.34)  for  a  and  p  to  solve  for  the  three  parameters  of 
the  bubble  motion:   rQ,  K  and  y« 

Suppose,  for  example,  that  we  know  P  when  the  momen- 
tiim,  s,  is  zero.   This  means  that  e  =  0  and  u  =  1.  Con- 
sider the  ratio  P/a.   We  have 


and  the  right  side  depends  only  on  k  and  y«   Combining 

o 

this  with  the  formula  for  a  /p ,  which  also  depends  only 
on  k  and  Yj  we  will  be  able  to  solve  for  these  quantities 
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Section  IV 
The  Electrostatic  Problem 


In  Section  III  functions  ^-,  ,  ^Cg  and  9?„  were  Intro- 
duced Into  the  energy  equation  to  take  account  of  the 
presence  of  free  and  rigid  surfaces.  However,  the  re- 
sults obtained  depended  only  on  the  ftmctlon  9?,  .   This 
function  9?,  can  be  evaluated,  to  a  first  approximation, 
by  finding  the  kinetic  energy  of  water  motion  due  to  a 
fixed  expanding  bubble.   Consider  the  velocity  potential 
for  a  fixed  expanding  bubble,  assumed  spherical,  located 
at  a  distance  H  beneath  the  s;arface  In  water  whose  depth 
Is  D  feet.   If  we  take  spherical  coordinates,  R,  9,  ^ 


H 


® 


with  pole  at  the  center  of  the  bubble,  the  problem  can  be 
formulated  mathematically  as  follows : 

Find  a  potential  function  f (R,e)  such  that 


(4.1) 
(4.2) 


•^   =  A   on  the  sphere, 
on. 

d$ 

•i-=-  =  0   on  bottom, 

on 


(4.3) 


$  =  0   on  surface. 
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The  last  condition  is  obtained  by  taking  Bernoulli's 
equation  on  the  surface,  assvming  that  gravity  can  be 
neglected,  and  then  neglecting  square  terras. 

Eqxiation  (4.1)  is  satisfied  by  taking 

This  is  equivalent  to  assuming  a  source  of  strength  A  A  at 
the  center  of  the  sphere.   We  can  then  satisfy  equations 
(4.2)  and  (4.3)  by  reflecting  this  source  in  the  bottom 
and  in  the  free  surface.   It  turns  out  that  the  image  with 
respect  to  one  boundary  must  be  reflected  in  the  other 
boundary  and  this  process  carried  out  on  the  successive 
images  leads  to  an  infinite  sequence  of  images:   Let 

(4.4)  $  =  A^A(^  +  $') 

where  $   is  the  potential  due  to  the  images.   It  can  be 
evaluated  by  a  method  similar  to  that  in  Report  37. IR. 

After  $  has  been  found  we  get  the  kinetic  energy  of 
the  water  by  evaluating 


|p/*i^ 


over  all  boundaries.   Because  of  conditions  (4.2)  and 
(4.3)  this  reduces  to  the  integral  over  the  sphere. 
Using  equation  (4.1)  we  have  that 

(4.5)   E=-gP    /   ^||dS=A  J$A^dco  =  A^A  /  fdco 

sphere 

2 
where  dco  is  the  element  of  solid  angle  so  that  ds  =  A  ^<*^- 

By  the  Mean  Value  theorem 

^  J$'  dco=  $'(0) 
where  $  (0)  is  the  value  of  the  potential  $  at  the  center 
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of  the  sphere.  Prom  equation  (4.4)  and  this  result  we 
rind  that 


(4.6) 


E  =  2TrpA^A^(|  +  $'(0)) 
=  SirpA^A^d  +  A$'(0)) 


so  that  <p-^   =  ^    (0) . 

Solution  of  the  Electrostatic  Problem 

This  shows  that  we  need  the  value  of  $  only  at  the 
center  of  the  sphere.   It  will  be  found  by  considering 
the  infinite  set  of  images  produced  by  reflecting  the 
original  source  at  the  center  and  its  images  in  the  bot- 
tom and  free  surface. 

For  convenience  of  representation  we  shall  draw  the 
svirfaces  vertically.   The  situation  is  as  follows: 


'-4 


^-3  ^-2 


'-1 


o 

S 

o 

O 

a 

-p 

^ 

-p 
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O 

m 
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® 

^ 

Oh 

■ 

© 

^O 

31^2 


Sg     S4 


Here,  S   is  the  source  of  strength  one  at  the  center 


Of  the  bubble.   S-j_,  S_-^,    Sg,  S_2,  •••  denote  the  image 
soTJirces  and  sinks.   They  are  obtained  as  follows: 
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S-,    is  the  reflection  of  S^   in  the  bottom 
S  ^   is  the  reflection  of  S    in  the  surface 


b 


-,  is  the  reflection  of  S  _   in  the  botton 
n+1  -n 


S    is  the  reflection  of  S   n  in  the  surface, 
-n  n-i 

llote  that  a  reflection  in  the  bottom,  to  satisfy  equation 
(4.2),  gives  an  image  of  the  saine  strength,  while  a  reflection 
in  the  free  surface  to  satisfy  equation  (4.o)  gives  an  image 
of  the  same  strength  v;ith  the  opposite  sign.   (A  source  of 
negative  strength  is  to  be  considered  a  sink.)   V.'e  find  that 

tiie  strength  of  Sp     is  (-)^   at  a  distance  2nD    from  S  ■ 

the  strength  of  Spj^.-i   is  (-)"   at  a  distance  2nD-2B  from  S 

the  strength  of  3  o    is  (-)    at  a  distance  2nD    from  S^ 
-^       -2n  o 

the  strength  °^  S  p  _i  ^^  (-)    at  a  distance  2nD+2B  from  S  . 

$  (0),  the  potential  at  the  center  of  the  bubble  due  to 
this  collection  of  sources  and  sinks  can  be  obtained  by  com- 
bining the  effects  of  S  and  S   . 

The  potential  at  S   due  to  Sp  and  S_p  is 


nD 


v/hile   that  due   to   S.^      •,    and  S_o     -,    is 

(4  7)       (_)^i L L_T  = (-)^    (43-2D) 

^    '    '       ^    '    '■(2n+2)D-2B     2nD+2B^       [  (2n+l)D+D-2B]  [  (2n+l)D-(D-2B)  ] 

The  total  potential  due  to  all  the  image  sources  and  sinks 
is,  therefore, 

00  /  xn   -,  .„      00  /  \n 

(4.3)     Z  i^  +  i(f  -  2)  Z: J^ g^-g 

1  ^    D  ^       0  (2n  +  1)2  -  (1  -  •^)2 
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2R 

or,    if  we   let  ^   -   1  =  x,    this   can  be  v/rltten  as 

(4.9)  .  iSp  .  f  g L-J.; . 

^  ^     0     (2n  +  1)"^   -  tl'^ 

The  series  on  the  right  can  be  expressed  in  teims  of 
tabulated  functions  as  follows: 

f  (X)  =  2  x: ^"^  / P 

0  (2n  +  1)'^  -  x'^ 
(4.10) 

=  2  ^[log  n— ^~)  +  /  (— ^-)]  -  ^  tan  -^. 

and  we  have 

$'(0)  =  -  i(f(x)  +  log  2) 
(4.11) 

Note,  that  if  D  — >  oo ,  $' (0)  =  -  -gg  .   By  differentiation 
we  find  that 

I 

/  /t  T  o  \        9$   _    1  ^^ '  /  „  \  9x  _     2   ~ '  ,  V 
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Appendix 


Table   3 


Table   4 


X 

f(x) 

0 

0 

.05 

.09184 

.10 

.1852 

.15 

.2816 

.20 

.3829 

.25 

.4910 

.30 

.6081 

.35 

.7579 

.40 

.8840 

.45 

1.051 

.50 

1.246 

.55 

1.481 

.60 

1.769 

.65 

2.156 

.70 

2.619 

.75 

3.292 

.80 

4.298 

.85 

5.968 

.90 

9.305 

.95 

19.31 

1.00 

00 

X 

^(x) 

X 

^(x) 

0 

.3466 

.05 

.3729 

-.05 

.  3157 

.10 

.3952 

-.10 

.2794 

.15 

.4143 

-.15 

.2366 

.20 

.4304 

-.20 

.1862 

.25 

.4440 

-.25 

.1264 

.30 

.4554 

-.30 

.05528 

.35 

.4651 

-.35 

-.03021 

.40 

.4751 

-.40 

-.1336 

.45 

.4796 

-.45 

-.2593 

.50 

.4849 

-.50 

-.4150 

.55 

.4891 

-.55 

-.6104 

.60 

.4924 

-.60 

-.8608 

.65 

.4950 

-.65 

-1.190 

.70 

.4969 

-.70 

-1.637 

.75 

.4982 

-.75 

-2.274 

.80 

.4991 

-.80 

-3.244 

.85 

.4996 

-.85 

-4.880 

.90 

.4999 

-.90 

-8.181 

.95 

.49998 

-.95 

-18.15 

1.00 

.5 

-1.00 

-00 
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